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ABSTRACT 

We analyze electron flux maps based on RHESSI hard X-ray imaging spec- 
troscopy data for a number of extended coronal loop flare events. For each 
event, we determine the variation of the characteristic loop length L with elec- 
tron energy E, and we fit this observed behavior with models that incorporate 
an extended acceleration region and an exterior "propagation" region, and which 
may include collisional modification of the accelerated electron spectrum inside 
the acceleration region. The models are characterized by two parameters: the 
plasma density n in, and the longitudinal extent Lq of, the acceleration region. 
Determination of the best-fit values of these parameters permits inference of the 
volume that encompasses the acceleration region and of the total number of par- 
ticles within it. It is then straightforward to compute values for the emission 
filling factor and for the specific acceleration rate (electrons s _1 per ambient elec- 
tron above a chosen reference energy). For the 24 events studied, the range of 
inferred filling factors is consistent with a value of unity. The inferred mean 
value of the specific acceleration rate above Eq = 20 keV is ~10~ 2 s" 1 , with a la 
spread of about a half-order-of-magnitude above and below this value. We com- 
pare these values with the predictions of several models, including acceleration 
by large-scale, weak (sub-Dreicer) fields, by strong (super-Dreicer) electric fields 
in a reconnecting current sheet, and by stochastic acceleration processes. 

Subject headings: Acceleration of particles — Sun: flares — Sun: X-rays and 
gamma-rays 
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Introduction 



The Ramaty High Energ y Solar Spectro s copic Imager (RHESSI, 



revealed a new class of flares (ISui et al.l 12004 ; IVeronig fc Brownl 12004 



Lin et al. 


200i 


!) has 


Krucker et al. 


2008; 



Xu et al.l 120081 ) in which the bulk of the hard X-ray emission is produced predominantly 



throughout the coronal portion of th e loop, but for which the hard X-ray spectrum, and the 



variation of source size with energy (IXu et al.l 120081 ). indicate strongly that such sources are 



not thermal. Rather they involve the injection of nonthermal electrons into a coronal region 
which is not only the site of particle acceleration, but also dense enough to act as a thick 
target, stopping the accelerated electr ons before t hey ca n penetrate to the chromosphere. 



As discussed by Ku et al.l (120081 ) and iGuo et al.l ( l2012al ). analysis of the variation of the 



longitudinal (field-aligned) source extent L of such sources, as a function either of photon 
energy e (deduced from im ages of the hard X -ray flux) or of electron energy E (deduced 



from electron flux maps - iPiana et al.l 120071 ) can be used to obtain information on the 



size of the acceleration region and the dens ity in the region into which the accelerated 



electrons propagate. Further, as discussed in IGuo et al.l ( I2012bl ). this information, coupled 



with information obtained through parametric fitting of the spatially-integrated hard X-ray 
spectrum, can be used to obtain information on the filling factor (the ratio of hard-X- 
ray-emitting volume to the total volume encompassing the source (s)) and on the specific 
acceleration rate, the ratio of the rate (s -1 ) of acceleration of electrons above a chosen 
reference energy to the number of electrons available for acceleration. 

While it is clear that collisional modification to the accelerated electrons occurs in the 
"propagation" region exterior to the acceleration region, an interesting question is: "should 
collisional modification to the accelerated electron spectrum inside the acceleration region 
itself also be considered?" On the one hand, the presence of a finite ambient density in the 
acceleration region should require inclusion of such a term. On the other hand, it could 
be argued either (1) that the acceleration mechanism "swamps" any effects of Coulomb 
collisions, or (2) that collisional processes are already included in the physics governing 
the form of the accelerated electron spectrum, so that an additional "post-acceleration" 
collisional term need not be considered. In either of these latter cases, we should model the 
acceleration region by an effectively "tenuous" region which has zero density for the purposes 
of calculating post-acceleration collisional modifications to the electron spectrum, although it 
must be stressed that the actual density in the acceleration region is likely to be comparable 
to that in the exterior "propagation" region, so that substantial bremsstrahlung emission 
will still be produced in the acceleration region. 

Alth ough "tenuous" acceleration region models have been considered by a n umber of 
authors jEmslie et al.l 120081 : IXu et al.l I2OO8I : iKontar et all boilL IGuo et al.l l2012al Jbh , so far 
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only lXu et al 



(200 8[) have consid ered the "dense" acceleration regiom model. Here we extend 



Guo et al.l ( 12012bl ) to include both "tenuous" and "dense" acceleration region 



the results of 
scenarios. 

In Section [2] we present the list of 24 extended-loop events studied. In Section [3] we 
derive analytical forms for L(E), both for the "tenuous" acceleration region model and for 
two different "dense" acceleration region models, one with a spatially uniform injection term 
over a finite length and one with a Gaussian spatial profile of the injection term. Fitting all 
three of these forms to the observationally-inferred forms of L(E) results in corresponding 
best-fit estimates of the acceleration region length Lq and density n. These are then used 
(Section H]) to infer, for each event, the value of the emission filling factor / (the fraction of 
the observed volume in which hard X-ray emission occurs). This is in turn used to estimate 
the acceleration region filling factor / acc (the ratio of the volume that is actively involved 
in electron acceleration to the overall volume that encompasses the acceleration region(s)), 
and so to obtain values of the specific acceleration rate 7]{> Eq), the rate of acceleration 
of electrons to energies > Eq divided by the number of ambient electrons available for 
acceleration. In Section we summarize the salient features of the results obtained and 
in Section O we compare the ensemble- averaged values of the quantity rj (> Eq) with the 
predictions of various acceleration models. In Section [7| we summarize the results obtained 
and present our conclusions. 



2. Observations 



The l is t of ev ent studied is shown in Table [TJ these are the same as those used by 
Guo et al.l (j2012bl ). with the addition of two additional events associated with a flare on 
2005 August 23 (Event #s 23 and 24). For each event, we fit the spatially-integrated hard 
X-ray emission spectrum with an isothermal-plus-power-law form, yielding values of the 
emission measure EM (cm -3 ) and temperature T (keV) of the thermal source, the intensity 
and spectral index 5 = 7 + 1 of the injected nonthermal electron spectrum (corresponding to 
the hard X-ray spectral index 7), and E t (keV), the transition energy between the thermal 
and nonth ermal compon ents of the hard X-ray spectrum. Straightforward thick-target mod- 



eling (e.g., lBrownlll97ll ) then provides values of dj\f/dt(> Eq) (s 



of electrons to energies above the (arbitrary) reference energy Eq 



the rate of acceleration 
20 keV. 



5 In this context, an "event" is a time interval during a flare for which spatial and spectral observations 
are sufficiently good to permit both a determination of the source spatial structure at a variety of energies 
and the overall spectrum of the hard X-ray emission. As can be seen in Table [TJ some flares provide multiple 
"events" ; other flares only one. 
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Table 1. Event List and Spectral Fit Parameters 


Event No. 


Date 


Time (UT) 


EM (10 49 cm" 3 ) 


T (keV) 


6 


E t (keV) 


dM/dt (10 35 s- 1 ) 


1 


2002-04-12 


17:42:00-17:44:32 


0.30 


1.53 


8.24 


15.5 


2.71 


2 




17:45:32-17:48:00 


0.46 


1.54 


8.01 


15.5 


4.69 


3 


2002-04-15 


00:00:00-00:05:00 


0.22 


1.75 


7.48 


15.5 


4.70 


4 




00:05:00-00:10:00 


0.76 


1.61 


7.93 


15.5 


9.32 


5 




00:10:00-00:15:00 


1.02 


1.60 


8.37 


15.5 


11.41 


6 


2002-04-17 


16:54:00-16:56:00 


0.06 


1.51 


5.70 


15.5 


0.39 


7 




16:56:00-16:58:00 


0.22 


1.43 


8.78 


14.8 


2.43 


8 


2003-06-17 


22:46:00-22:48:00 


1.92 


1.71 


9.95 


16.5 


17.27 


9 




22:48:00-22:50:00 


2.59 


1.67 


10.36 


16.5 


17.91 


10 


2003-07-10 


14:14:00-14:16:00 


1.26 


1.45 


10.05 


15.5 


7.43 


11 




14:16:00-14:18:00 


1.31 


1.34 


10.38 


14.8 


8.53 


12 


2004-05-21 


23:47:00-23:50:00 


0.35 


1.85 


7.07 


18.5 


3.28 


13 




23:50:00-23:53:00 


0.62 


1.75 


7.51 


18.5 


2.32 


14 


2004-08-31 


05:31:00-05:33:00 


0.06 


1.61 


10.56 


15.5 


0.40 


15 




05:33:00-05:35:00 


0.21 


1.57 


12.19 


18.5 


0.29 


16 




05:35:00-05:37:00 


0.29 


1.48 


7.45 


18.5 


0.16 


17 


2005-06-01 


02:40:20-02:42:00 


0.14 


1.81 


6.53 


17.5 


1.44 


18 




02:42:00-02:44:00 


0.37 


1.70 


7.86 


17.5 


2.67 


19 


2011-02-13 


17:33:00-17:34:00 


0.54 


1.39 


5.86 


10.5 


35.02 


20 




17:34:00-17:35:00 


0.52 


1.68 


6.55 


14.5 


19.43 


21 


2011-08-03 


04:31:12-04:33:00 


0.36 


1.61 


9.23 


15.5 


3.96 


22 


2011-09-25 


03:30:36-03:32:00 


0.13 


1.44 


8.33 


14.5 


1.19 


23 


2005-08-23 


14:23:00-14:27:00 


0.04 


2.12 


6.32 


14.2 


5.44 


24 




14:27:00-14:31:00 


0.29 


1.87 


7.96 


16.5 


10.20 



-5 - 



For each event, we also created electron flux images at a series of e lectron energies E , 
each produced using electron visibilities constructed via the procedure of IPiana et al. ( 2007 ) 
and the uv_smooth image reconstruction algorithm ( Massone et al. 20091 ). At each energy 
E, the (field-aligned) root-mean-square len gth L and (cross-fi eld) width W of the source 
were then calculated using the procedure of iGuo et al.l (12012al ). 



3. Analysis 

3.1. Tenuous Extended Acceleration Region Model 



model flXu et al. 



Guo et al.l ( 2012b!) fit the observed form of L(E) to a "tenuous" acceleration region 



2008 



j). In such a model, either (a) the acceleration process is considered to 



dominate the evolution of the electron spectrum to such an extent that the effect of collisions 
in the acceleration region can be ignored, or (b) the effect of collisions is already incorpo- 
rated in the production of the accelerated electron spectrum. In either case, no additional 
collisional modification to the electron spectrum need be considered. Thus, the "tenuous" 
acceleration region is considered to have zero density only in respect of the variation of the 
electron spectrum throughout the acceleration region. However, it must again be stressed 
that in such a model there is still a finite density (and so a bremsstrahlung target) in the 
acceleration region, so that the hard X-ray source includes the acceleration region itself. 

Consider, then, electrons with an injected spectrum F (E ) ~ Eq 5 (cm -2 s -1 keV" 1 ) 
that are accelerated within a region extending over [— Lo,t/2,Lo,t/2] an d propagate through 
an exterior region with uniform density n t , in which the electrons suffer energy loss through 



Coulomb collisions with ambient particles, particularly electrons ( lEmslidll978l ). Since col- 
lisions in the acceleration region are neglected, the calculated value L t (E) in such a model 
consists of two parts: the (energy-independent) length of the acceleration region L o t and 
the (injection-spectrum-weighted) average penetration depth associated with the electrons 
of injection energies E$> E that contribute to the electron map at energy E. For Coulomb 
collisions in a uniform ambient target, the latter term is equal to a spectral shape factor 
times the quantity E 2 /2Kn t , where K = 27re 4 A, e being the electronic charge and A the 



Coulom b logarithm. Quantitatively, L t (E) has the simple quadratic form (see IGuo et al. 



2012bh : 



Least-squares fitting of the observed values of L(E) to this analytic form leads to best-fit 
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values of the acceleration region length Lo.t and the density n t in the exterior propagation 
region, which we also take to be the dens ity in the accelerat ion region itself. The resulting 
values of L o t and n t have been given in iGuo et al.l (|2012bl ). and, for completeness, these 
results are repeated here in Table 2. 



3.2. Dense Extended Acceleration Region Model 

As noted above, the form (Tj[|) does not include the effects of collisional modification 
to the accelerated electron spe ctrum within th e acceleration region itself. Therefore, for 



completeness we here, following IXu et al.l ( 120081 ). develop the expression for L(E) associated 



with a "dense" acceleration region model, which does take into account post-acceleration 
collisional losses within the acceleration region itself. For the tenuous acceleration model, the 
lack of collisional modification to the electrons inside the acceleration region means that the 
form for L t (E) - Equation ([I]) - does not depend on the form of the variation of the electron 
injection rate with position in the acceleration region. However, for the dense acceleration 
model, collisional modification to the electron spectrum in the acceleration region means 
that the results do depend on the spatial form of the acceleration rate. Hence, we have 
evaluated the form of L(E) for two illustrative parametric forms of the spatial profile of the 
accelerated electron flux, in order to evaluate the sensitivity of the results to the particular 
form of this spatial profile. 



3.2.1. Spatially Uniform Injection Profile 

We first consider an acceleration region extending over [— Lo )W /2,Lo )lt /2], with a uniform 
density n u equal to that in the rest of the loop and standard collisional losses appropriate 
to such a density. Electrons, with a spectrum F (E ) ~ Eq S , are assumed to be injected 
uniformly throughout this region, and the electrons injected at each location s (— L , u /2 < 
s o < Lo jU /2) propagate in both directions along the loop. The electron flux spectrum at 
longitudinal position s then follows from the one-dimensional continuity equation: 



F U (E, s, s ) = F (E ) — ~ {E2 + 2Knuls _ Sol){5+m » ( 2 ) 

where we have used the collisional relation Eq = E 2 + 2KN = E 2 + 2Kn u \s — s a \, N being 
the column density (cm -2 ) along the direction of electron propagation. The corresponding 
expression for the longitudinal loop extent L(E) is: 
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L U (E) 



(3) 



Note that we have formally approximated the limits on the s-integration at ±00, rather than 
using an estimated overall loop extent ±L max . This approximation was used principally to 
avoid having to estimate the total loop length, including the portion occupied by weakly- 
emitting regions near the chromosphere which, especially for events far from the limb, are 
the regions where the observed geometry is most sensitive to projection effects. However, 
since L U (E) is equal to the square root of a ratio of integrals, the results do not depend 
significantly on the values of the s-integration limits used. 

The forms ([2]) and Q do not result in a closed analytic form for L U (E) (cf. Equation ([ID) 
and so a straightforward least-squares fitting procedure cannot be used to determine the best- 
fit values of the acceleration region length Lo,« and density n u . However, these best-fit values 
can nevertheless be determined numerically by exploring (L 0jU ,n u ) space and computing, at 
each sampled point, the value of x«(= Yl^T ax W (E)[L U (E) — L(E)} 2 /rif ree ), where the point 
weighting^ w(E) = l/L 2 err {E) and nj ree is the number of degrees of freedom for the fitting. 
The best-fit values of Lq )U and n u are determined by finding the location of the minimum 
value of xi- 



The left panel of Figure [T] shows the Xu(^o,u,n u ) surface for Event #1 (2002 April 12), 
together with contours of xt equal to 1, 2, 5, and 10. A pronounced local minimum (xl, — 
0.13) at (L 0jU = 33 arcseconds, n u = 1.2 x 10 11 cm" 3 ) is found. By contrast, the right panel 
of Figured] shows the Xu(Lo, n u) surface for Event #9 (2003 June 17), together with contours 
of xt equal to 3, 5, and 10. For this event, the local minimum of and so the best-fit 
values (L 0)U = 35 arcseconds, n u — 1.0 x 10 11 cm -3 ) are, especially for the density n u , not 
so precisely determined, and this is reflected in the associated higher value xt = 2-07 in 
Table 2. 

The best-fit values of of L 0tU and n u , and the corresponding values of x\> f° r ai l 24 
events are listed in Table 2. 



In order to provide quantitative uncertainties L err (E) on observed loop lengths L(E), we applied a 
Monte Carlo approach in wh ich random noise is added to the visibilities and the resulting uv-smooth images 
recomputed and reanalyzed (|Guo et al.ll2012al) . 
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Fig. 1.— Left: The xl(L ,u,n u ) surface for Event #1 (2002 April 12). Right: The 
xi(LQ tU ,n u ) surface for Event #9 (2003 June 17). The minimum values of x\ are a l so 
shown. 



3.2.2. Gaussian Injection Profile 

Here we again assume an extended, dense acceleration region, but for which the spatial 
distribution of the energetic electron injection falls off with distance sq from the "kernel" of 
the acceleration site (of uniform density %) according to a Gaussian form characterized by 
a standard deviation cr . For such an injection profile, the form of the electron flux spectrum 
F G (E, s, s ) is 

F G (E s sq) ~ ex P(- s o 2 / 2ff o 2 ) E (4) 

G{ ' ' 0j (T (E 2 + 2Kn G \s - So |) (m)/2 ' 1 ' 

which is substituted into 



L G (E) = f-oof-oo s2FG ( E i s > s ^ ds ° ds (5 ) 
° V lZof-oo F G( E ^^o)dsods 

to determine the form of the longitudinal loop extent L G (E). The best-fit values for o~ Q and 
n G can then be derived by a procedure similar to that for the uniform injection case. 



The left panel of Figure [2] shows the Xci^o^c) surface for Event #1 (2002 April 12), 
together with contours of x G equal to 1, 2, 5 and 10. A pronounced local minimum (x G = 
0.12) at (ctq = 13.5 arcseconds, n G = 1.2 x 10 11 cm -3 ) is found. By contrast, the right 



Fig. 2.— Left: The XgOo,™g) surface obtained for Event #1 (2002 April 12). Right: The 
Xg^o^g) surface for Event #9 (2003 June 17). The minimum values of \ G are also shown. 

panel of Figure |2] shows the Xg^o^g) surface for Event #9 (2003 June 17). Again, for the 
latter event, the local minimum of Xg ( an d so the best-fit values of a = 17.5 arcseconds 
and ric = 0.9 x 10 11 cm -3 ) are not so precisely determined. This is especially true for the 
value of the characteristic scale ctq, and this lack of precision is reflected in the associated 
high value Xg = 2-05 in Table 2. 

We define the characteristic length of the acceleration region as the full width of the 
Gaussian function at half the maximum (FWHM) L , G = 2 V2 \n2a . The best-fit results 
for Lo,g and Ug, and the corresponding values of Xg-> f° r an 24 events are shown in Table 2. 

In general, the values of both L and n for the Gaussian injection profile are compa- 
rable to those for the uniform injection case, illustrating that the results are not critically 
dependent on the spatial profile of the acceleration rate within the acceleration region. 

4. Derived Quantities 

4.1. Emission Filling Factor 

The soft X-ray emission measure EM is related to the plasma density n and the emitting 
volume Kmit through EM = n 2 V cm \ t - Given that an emitting region (considered here to be 
a cylinder of length L and width W) may be composed of a number of discrete emitting 
subregions (e.g., "strands," "kernels"), the emitting volume may be equal to or smaller than 
the total flare volume 
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V = -W 2 L(E t ), (6) 

a quantity readily ascertained from RHESSI images at the transition energy E t (Tabled]), 
the maximum energy at which thermal emission is predominant. 

The ratio of the emitting volume Kmit, determined from the values of the emission 
measure EM (Table [1]) and the inferred acceleration region density (n u or no, as appropriate; 
see Table 2), to the volume V of the observed region that encompasses the emitting region(s) 
is termed the emission filling factor (or simply "filling factor"): 

- = ]4mt = EM 

; V n 2 V ' 1 ' 

Values of / for each event, and the geometric mean and multiplicative la uncertainty for 
the ensemble of events studied, are given in Table 2. 



4.2. Specific Acceleration Rate 

The specific a cceleration rate rj(> Eq) (electrons s _1 per ambient electron) is defined 



(jEmslie et al.ll2008l ) as the ratio of two quantities: dN /dt(> Eq), the rate of acceleration of 
electrons to energies greater than £"0, and A/", the number of particles available for accelera- 
tion: 

The quantity dAf/dt(>20 keV) is readily determined by spectral fitting of the spatially- 
integrated hard X-ray emission; values for all 24 events are given in Table [TJ The quantity 
M can be estimated from the relation 

ttW 2 

A/"= /acc^j— ^0™, (9) 

where the values of the acceleration region length Lq and density n are inferred from the 
appropriate acceleration region model and f acc is the acceleration region filling factor^ the 
ratio of the volume in which acceleration actually occurs to the encompassing volume of the 



7 Note that the ap pearance o f the fa ce term in the expression for Af was implicitly assumed to be unity 
in the car lier work of lGuo et all |2012bl ). Here we elaborate on this assumption further. 
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region that contains the acceleration site(s). This quantity is not to be confused with the 
emission filling factor /, since the regions in which acceleration occurs are not necessarily 
congruent to the regions that produce hard X-ray emission. 

We have no direct knowledge of the value of the acceleration region filling factor / acc . 
However, for the coronal-loop sources considered here, it appears reasonable to assume that 
the volume in which acceleration occurs is a strict subset of the acceleration-plus-propagation 
volume. Since hard X-ray emission occurs in both the acceleration and propagation regions, 
it follows that the acceleration filling factor is less than the emission filling factor, i.e., that 
face < f ■ Further, a value f acc that is less than unity reduces the number of electrons Af 
available for acceleration, and therefore, for a given acceleration rate dj\f/dt(> Eq), the 
specific acceleration rate r)(> E ) is increased. As we shall see, the values of rj(> E ) 
obtained for / acc = 1 already impose significant constraints on particle acceleration models, 
so that a value / acc <C 1 is extremely unlikely. 

Table 2 shows that the values of the emission filling factor / do vary considerably from 
event to event, and indeed that the values of / obtained for some events are greater than 
one (and thus unphysical). Such results are doubtless artifacts of vagarie^fl in the inferred 
electron maps and/or in the modeling assumptions used. Apart from a few unusual cases 
(discussed in Section [5]), the value of the emission filling factor / is generally bounded above 
by ~(2-3), well within the bounds of possible uncertainty in the quantitative electron flux 
values inferred from the RHESSI data and in the modeling assumptions used. Moreover, the 
geometric mean of the emission filling factor / for the ensemble of 24 events studied is very 
close to unity, a value hitherto tacitly assumed by many authors. We therefore assume that 
/ = 1 (i.e., that the entire observed region V participates in both the hard X-ray emission 
process) and, based on the arguments in the previous paragraph, we further make the rather 
bold assumption that the emission and acceleration regions are congruent, so that / acc = 1. 
In doing so, we recognize that this may be an overly constraining assumption for models 
in which the acceleration region is intrinsically a small fraction of the region in which the 
electrons propagate and emit bremsstrahlung (see, e.g., Section l£T2l) . 

From Equations (jSJ) and (jUJ), with f acc = 1, we can deduce the value of the specific 
acceleration rate rj(E ) in each event; values (for E = 20 keV), together with the geometric 
mean for the ensemble and its la multiplicative uncertainty, are given in Table 2. 



For example, the derivation of quantitative electron flux maps depends (jPiana et al on the form of 

the bremsstrahlung cross-secti on used. This cross-s ection can vary a great deal with the degree of directivity 
of the exciting electron beam ([Massone et al.ll2004l ). 
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Table 2. Acceleration Region Characteristics 



Event No. 




Lo,u 


Lq,g 


n t 


n u 


n G 


ft 


fu 


fc 


Vt 


l/u 


VG 


X 2 t 


X U 


X 2 a 






(arcsec) 




(10 


11 cm" 


~ 3 ) 








(> 20 keV; xlO" 


' 3 s" 1 ) 










J-O.U 


oo.u 


^1 8 
ol.o 


1 ^ 


1 2 


1 2 




u.uy 


u.uy 


6.5 


4.6 




12 


n i 1 ) 

U. -LD 


12 


2 


16.3 


31.0 


24.7 


1.4 


1.4 


1.3 


0.83 


0.82 


1.08 


14.5 


7.3 


10.5 


0.33 


0.37 


0.33 


3 


16.7 


29.0 


24.7 


4.4 


1.1 


1.3 


0.04 


0.56 


0.47 


4.0 


8.9 


9.5 


0.31 


0.34 


0.26 


4 


16.6 


31.0 


24.7 


4.8 


1.7 


1.5 


0.11 


0.93 


1.12 


7.3 


11.4 


15.6 


0.34 


0.41 


0.28 


5 


16.6 


29.0 


22.4 


10.5 


2.6 


2.3 


0.03 


0.43 


0.56 


3.3 


7.5 


11.1 


0.41 


0.52 


0.82 


6 


11.9 


19.0 


13.0 


4.9 


1.0 


0.9 


0.02 


0.39 


0.43 


0.6 


2.1 


3.2 


0.14 


0.34 


0.36 


7 


10.4 


19.0 


15.3 


1.8 


0.8 


0.9 


0.44 


2.41 


1.83 


12.1 


15.5 


16.7 


0.18 


0.24 


0.26 


8 


17.8 


33.0 


29.4 


2.6 


0.8 


0.9 


0.90 


8.99 


8.20 


24.1 
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5. Results 

Results for all events are shown in Table 2. The main salient features are as follows: 

• The loop densities n u and na from each of the dense acceleration region models are 
typically smaller, by a factor ~3, than the value n t deduced for the tenuous acceleration 
region model. This is because in dense models collisional energy loss also occurs within 
the acceleration region itself, so that the overall density does not have to be as high 
to accomplish attenuation of the electron flux over the overall distance that is evident 
from the electron flux images; 

• The sizes of the acceleration region Lo iU and Lo,g m dense acceleration region models 
are generally somewhat larger than the acceleration region lengths L Qtt deduced from 
the tenuous acceleration region model. 

• The filling factors f u and fa for the dense acceleration region models are about an 
order of magnitude larger than the values of ft obtained from the tenuous model. 
Further, both and fa are consistent with (the upper limit value of) unity, albeit 
with a large (lex) spread of a factor of about four. For some events (e.g., Events 8, 9, 
10, 11 and 15), the values of both f u and fa are in excess of unity by an unacceptably 
large factor. These events are also characterized by unreasonably large values of the 
acceleration region length Lq (comparable to, or even larger than, the overall extent 
of the observed loop) and, in general, the x 2 values of these fittings are unacceptably 
large (see Table 2), so that these values of / and Lq should not be taken too seriously. 

• The values of the specific acceleration rate ry(>20 keV) for both of the "dense" acceler- 
ation region models are about a factor of two larger than in the "tenuous" acceleration 
region model. This is because the neglect of collisional losses in the acceleration region 
(effectively modeled by assuming a lower density within the acceleration region) im- 
plies a smaller number of total particles M in the acceleration region, so that a given 
acceleration rate dN /dt(> E ) (deduced from observations of the spatially-integrated 
hard X-ray flux) corresponds to a larger specific acceleration rate. The geometric mean 
values for ?7(>20 keV) are ~5xl0 -3 s _1 for the tenuous acceleration region model and 
~10xl0 -3 s _1 for the dense acceleration region model (for both uniform and Gaussian 
spatial injection profiles). Each of these values has a la spread of about a factor of 
three. 

• The x 2 goodness-of-fit measures for all three models are comparable. This indicates 
that not only is the detailed spatial profile of the electron acceleration poorly con- 
strained by the observations, but also that the decision whether to include collisional 
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modification to the electron spectrum in the acceleration region itself is not driven by 
the data, but rather by the theoretical considerations related to the physical consis- 
tency of the overall scenario. Given this, it is also possible that a "hybrid" model, in 
which there is a difference (or longitudinal gradient) in the collisional energy loss term 
(effectively modeled by a spatially varying density) is also a viable model. 



Note that some o f the flares studied are located at significant distances from the solar 
limb (see Figure 2 in iGuo et all l2012bl ) and consequently the inferred loop lengths may 
be underestimated due to line-of-sight foreshortening effects. Further, given that the loop 
lengths generally increase with energy, the apparent source length L, and so volume V, 
at higher energies generally involves a greater portion of the foreshortened parts of the 
magnetic loop, and thus the foreshortening effect may be more significant at high energies 
than at lower energies. According to Equations through Q, underestimating the loop 
lengths L (and hence volumes V) will tend to result in overestimation of both the filling 
factor / and the specific acceleration rate rj. However, foreshortening effects also reduce the 
coefficient of E 2 in Equation (TO , causing the inferred loop densities n t to be overestimated 
in the tenuous acceleration model; similar results hold for the expressions for Ljj and Lq in 
the dense acceleration model - Equations ([3]) and According to Equations ([7]) and 
overestimating the density n tends to produce the opposite effect of an increased volume V, 
i.e., smaller values of both the filling factor / and the specific acceleration rate rj. Thus, 
overall, these two effects offset, so that the inferred values of the filling factor and specific 
acceleration rate are relatively insensitive to the effects of foreshortening. Nevertheless, it 
must be acknowledged that the effects of foreshortening are a possible source of uncertainty 
in our results. 



6. Implications for Acceleration Models 



In this section we discuss the significance of the results obtained, with particular atten- 
tion to the values of the specific acceleration rate 77. Specifically, we compare the values with 
the predictions of candidate electron acceleration models, respectively iny olving acceleration 



by large-scale sub-Dreicer electric fields (see, e. g..lBenka &: Holmanl 1994T) , by super-Dreicer 



electric fields in a reconn ecting current sheet ( jLitvinenko fc Somovl Il993l ; iLitvinenkol 11996 
Litvinenko fc Craig 120001 ). and by stochastic acceleration processes. 
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Fig. 3. — The variation of the specific acceleration rate rj with electric field strength £ (in 
units of the Dreicer field £d). The solid and dashed lines correspond respectively to the 
geometric mean value rj ~ 1.0 x 10~ 2 (electrons s _1 per ambient electron) in the dense 
acceleration region model and its la ranges (x/-r 3.0) as shown in Table 2. 



6.1. Sub-Dreicer Acceleration 



Application of a weak, sub-Dreicer, electric field £ over an extended distance leads to 
runaway acceleration of those electrons in the high-energy tail of the background thermal 
distribution for which the applied electric force exceeds the initial Coulomb drag force. An 
e mpirical expression f or the specific acceleration r ate r? i n such a scen ario has been give n 
bv lEmslie et aD toO^ ) [after iKruskal fc Bernsteml Jl964h . ICohe^l Jl976h . andlSinghj Jl977h ]. 



viz. 



exp 



£d 
A£ 



2£b 
£ 



(10) 



where £d is the Dreicer field ~ 2 x 10~ 8 n (cm~ 3 ) / T (K) V cm -1 . For the typical densities 
n ~ 10 11 cm -3 and temperatures T ~2 keV ~ 2 x 10 7 K involved (see Tables [T] and 2), 
£d ~ 10~ 4 V cm -1 . As shown in Figure El the variation of rj with £ is very strong; rj varies 
by more than fifteen orders of magnitude over the range £/£n = 0.01 to 1. 

For the dense acceleration region model, the average value of rj for the 24 events studied 
(Table 2) is ~ 1.0 x 10~ 2 , with a multiplicative uncertainty of about 3. (As noted in 
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Section [6], these values do not depend critically on the spatial form of the injection profile.) 
These bounds on 77, and the corresponding bounds on £/£&, are shown in Figure [31 from 
which it is apparent that for the sub-Dreicer model to account for the observed range of 77 
values, the applied electric field £ must be in the very narrow range from about 0.15 £d to 
0.3 £dj with an average value ~ 0.2 Sd- Such an electric field strength £ ~ 2 x 10~ 5 V cm -1 
rather nicely accelerates electrons up to the threshold energy ~20 keV, but only by using all 
the available (10 — 15) arcsecon d (~ 1 9 cm) half-length of the coronal loop. Hence, contrary 



to claims by iBenka fc Holmanl ( 119941 ) . this mechanism cannot account for the appearance 
of a power-law spectrum of accelerated electrons up to much higher energies ^ 100 keV. 
Further, it is extremely unlikely that the electric field strength (in units of the Dreicer field) 
lies within a factor of two in all 24 events studied. 

A similarly narrow range of electric field values is obtained if a tenuous acceleration 
region model is assumed to hold for all events, albeit with a slightly smaller value of 8/ Ed- 
Indeed, even allowing for a variation between tenuous and dense acceleration region models 
from event to event, the allowable range of £/£d values is very small. We therefore re- 
ject the hypothesis that acceleration by large-scale sub-Dreicer fields is responsible for the 
acceleration of electrons to hard- X-ray-producing energies in the events studied. 



6.2. Super-Dreicer Acceleration 



In this model (ILitvinenko &: Somovlll993l ; iLitvinenkol Il996l ; iLitvinenko fc Craigl 120001 ). 
the electric field £ = (v/c) x B is generated inductively by advective motions in the vicinity 
of the current sheet associated with the reversal of the principal component of the mag- 
netic field. For field strengths B ~ 100 G and inflow speeds of order the Alfven speed 
va = Bj -»/ 4:7inm p ~ 7 x 10 7 cm s _1 (where m p = 1.67 x 10~ 24 g is the proton mass), £ is 
of order 100 V cm -1 , far in excess of the Dreicer field. Therefore collisional effects in the 
core acceleration region itself are negligible; after being accelerated for only a few meters, the 
electrons reach hard- X-ray-producing energies of order 100 keV, and the collisional modifica- 
tion to the accelerated electron spectrum that defines the "dense acceleration region model" 
actually occurs between current sheets in the overall volume occupied by such sheets. Before 
the accelerated electrons can reach even higher energies, they escape the current sheet be- 
cause of the presence of a perpendicular magnetic field component that directs the particles 
out o f the s heet. The sp ecific acceleration rate in such a model is given by (jEmslie et al. 
2008h [after ISomovi dl992h ]: 



-17- 



,~(^f^i, di: 



where B± and B\\ are the perpendicular and parallel components of the magnetic field and 
L s is the length of the current sheet, which we take to be comparable to the inferred length 
Lq of the acceleration region containing the multiple current sheets. Thus, setting L s ~ 
Lq ~ 20 arcsec ~ 1.5 x 10 9 cm (Table 2), we find that, numerically, 



l? ~°- 1 ^J 8 ■ (12) 

The observationally-inferred typical value r/(>20 keV)~ 10~ 2 s _1 is thus obtained for a 
field aspect ratio B±/B\\ ~ 0.1. This is not an unreasonable value; furthermore, unlike for the 
sub-Dreicer model (in which the value of rj was extremely sensitive to the value of the applied 
electric field £), the value of rj scales only linearly with this aspect ratio. Consequently, the 
one-order-of-magnitude range of rj values corresponds to a very reasonable one-order-of- 
magnitude range in magnetic field aspect ratio: 0.03 ^ B±/B\\ ^0.3. We therefore conclude 
that super-Dreicer acceleration in a current sheet (or, more likely, set of current sheets) is a 
viable model for explaining the inferred values of the specific acceleration rate rj. However, 
there are clearly challenges associated with incorporation of such a current sheet geometry 
into the overall form of the observed hard X-ray images. Further, the inference of unit 
emission filling factors / (Section 14. ip . and the assignment of a unit acceleration region filling 
factor / acc (Section 14. 2[) present considerable geometric challenges for such an acceleration 
model. Moreover, as noted in Section 14.21 assigning a value / acc <C 1 in Equation ([9]) would 
result in a much smaller value of M and so, for a given value dM/dt{> E ), a value of 
v(> Eq) that is too high to be accounted for by Equation (I12p for reasonable values of the 
magnetic field aspect ratio (B±/B\\). 



6.3. Stochastic Acceleration 



Stochastic acceleration models are characterized not through the action of a coherent 
electric field acting over an extended distance, but rather through multiple interactions of 
the particle with scattering c enters, which may tak e the form of, for example, regions of 
enhanced plasma turbulence (IPetrosian fc Liul 120041) or enhanced magnetic field strength 
( [Miller et al.lll996[ ). In general ( iBian et al.ll2012[ ). the acceleration efficiency of a stochastic 
acceleration model is determined by the momentum diffusion coefficient D p = tlF 2 , where 
F 2 is the mean square magnitude of the accelerating force and tl is the force correlation 



- 18 - 



time in a Lagrangian frame that comoves with the accelerated particle. 



Equation (147) of iBian et al.l ( 120121 ) gives the evolution of the momentum phase-space 
distribution function f(p,t) under the action of a three-dimensional momentum diffusion 
process characterized by a momentum diffusion coefficient D(p) = Dop a : 



f(p,t)~t^ 2 -^ 



exp 



-p 



2-a 



(13) 



(2 -a) 2 D t_ 

Integrating from reference momentum p to oo shows that the number of particles accelerated 
above reference momentum p within time t is 



PO J XQ 



OO 2 — OL 



AT(po,t)~ / p 2 f(p,t)dt~ t W-^ / e- x x {l+a V {2 ~ a Ux- x = ^ 



(2-a) 2 D t 



(14) 

The dependence of M on p arises from the appearance of the lower limit Xq is the integral. 
Therefore, the characteristic rate (s _1 ) at which particles are accelerated to momenta greater 
than p Q depends on the magnitude of the diffusion coefficient D according to 

^(2-a) 2 A>^ 2 = (2-a) 2 ^. (15) 

Po 

With r\ ~ Dp/p 2 = tlF 2 /p 2 and writing F 2 = p 2 /r 2 cc , where r acc is a characteristic accelera- 
tion time, we find that the specific acceleration rate is related to the correlation time tl and 
the acceleration time r acc through 

T)- 1 ^- (16) 

T 

acc 

Moreover, the total number of scatterings that an electron undergoes can be estimated as 
the ratio of the acceleration time to the time spent in a single scattering event: 

Ns ~ 1^1 ~ — L . (17) 

TL V Tacc 

With observationally-inferred values rj ~ 10~ 2 s _1 , it follows that, if stochastic acceleration 
is the responsible mechanism, 



(18) 
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The acceleration time can be eq uated roughly with the rise time of the hard X-ray flux, 
typically a few seconds (see, e.g., lAschwanden et al.lll995l ). Adopting a value r acc = 3 s, we 
find that r L ~ 0.1 s and N s ~ 30. 



For stochastic acceleration by cascading magnetohydrodynamic turbulence, iMiller et al. 
( 11996k their Figures 6, 7, 9, 10 and 12) find a volumetric electron acceleration rate ~(1.5 
4) xlO 8 cm -3 s" 1 above 20 keV, with the exact val ue dependen t on th e parameters and 



assumptions of the various models considered. In the IMiller et al.l (119961 ) model, the back- 
ground number density is n = 10 10 cm" 3 , so that 77 ~(1.5 - 4) xl0~ 2 s" 1 above 20 keV. 
This value of 77 is somewhat higher than the values inferred from our analysis, although 
it must be noted that the value of the acceleration region density in their model is an or- 
der of magnitude lower than that inferred from our observations of extended coronal-loop 
sources. We therefore encourage modeling of stochastic acceleration in regions of density 
n ~ 10 11 cm -3 to determine values of 77 for comparison with our empirically-inferred values. 
For deka-keV electrons with velocity ~ 10 10 cm s -1 , the Lagrangian correlation time of order 
0.1 s corresponds to a correlation length ~ 10 9 cm, comparable to the length of the observed 
source. 



7. Summary and Conclusions 

We have shown that the structure of extended coronal-loop sources can be well explained 
both by an extended "tenuous" acceleration region model, in which the effect of collisions 
in the acceleration region is ignored (they are assumed either to be overcome by, or to 
be included in, the acceleration process) and by a "dense" acceleration region model, in 
which post-acceleration collisional energy losses affect the electrons equally inside and outside 
the acceleration region. The compatibility of both models with the observations therefore 
clearly allows for the possibility of a "hybrid" model, in which the effect of collisions varies 
throughout the acceleration and propagation regions, either due to actual density gradients 
and/or to a spatially- varying dominance of the acceleration process over collisional losses. 

Adopting a "dense" acceleration region model with a homogeneous collisional loss (den- 
sity) profile both inside and outside the acceleration region results in the emission filling fac- 
tor for the hard-X-ray-emitting plasma being closer to unity than in the "tenuous" model. In 
both models, there are a few anomalous cases, in which unphysical values of / greater than 
unity are realized. Discarding these cases, and also adopting a value of unity for the accel- 
eration region filling factor, the average value of the specific acceleration rate (electrons s -1 
per ambient electron) is r/(>20 keV) ~ (5 — 10) x 10~ 3 s -1 . 
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Such a value of rj is difficult to account for in an acceleration models that invokes large- 
scale acceleration by a weak, sub-Dreicer, electric field. For such models to be valid, the 
ratio of the applied electric field £ to the Dreicer field 8-q must lie in a very close range, 
extending over less than a factor of two. Such an acceleration model is, however, consistent 
with the inferred filling factor / of order unity. 

The inferred values of ?7(>20 keV) are broadly consistent with the predictions of models 
that invoke large, super-Dreicer, fields formed inductively during the collapse of magnetic 
field lines in the vicinity of an abrupt field reversal, i.e., a current sheet. The value of 77 scales 
linearly with the "aspect ratio" of the magnetic field geometry, and it is not unreasonable 
to expect values of the sheet aspect ratio that extend over the required range from ~(0.03 - 
0.3). However, the issue of the overall geometry of a source composed of multiple thin current 
sheets remains a significant challenge to this model, particularly in view of the finding that 
the emission filling factor / — 1. 

The stochastic acceleration model is consistent with both an emission filling factor 
/ ~ 1 and the values obtained for the specific acceleration 77(>20 keV). How ever, given 



the c onsiderable number of degrees of freedom associated with such models (see iBian et al. 



20121 ). it must be admitted that achieving agreement with observationally-inferred values is 
not so much a verification of the model, but rather a way of setting bounds on the parameters 
that define it. We have, however, concluded that a viable stochastic model has a Lagrangian 
correlation time (distance) of order 0.1 s (10 9 cm). Further, during the several seconds of 
the acceleration process, the accelerated electron is involved in roughly 30 collisions with 
scattering centers. 

There exists very little literature on the predicted valu e of 77 for stochastic acceleration 



models. What little does exist (see, e.g.. iMiller et al.lll996l ) is based on parameters that are 



not comparable to the physical environment that characterizes the dense coronal loop sources 
discussed herein. We therefore encourage modelers to explore stochastic acceleration models 
in parameter regimes that are more aligned with the physical conditions in such sources. 
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